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Abstract The different timing results of the magnetar Swift J1822.3 — 1606 is analyzed 
and understood theoretically. It is pointed that different timing solutions are caused not 
only by timing noise, but also that the period derivative is decreasing after outburst. Both 
the decreasing period derivative and the large timing noise may be originated from wind 
braking of the magnetar. Future timing of Swift J1822.3 — 1606 will help us make clear 
whether its period derivative is decreasing with time or not. 
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1 INTRODUCTION 

Magnetars are peculiar pulsar-like objects. They are assumed to be neutron stars powered by strong 
magnetic field decay (Duncan & Thompson 1992). A neutron star is often confirmed as a magnetar if 
its surface dipole magnetic field is higher than the quantum critical field (£?qed = 4.4 x 10 13 G). The 
surface dipole magnetic field is calculated from the period and period derivative (assuming magnetic 
dipole braking, Kouveliotou et al. 1998). However, the magnetic dipole braking assumption will also 
provide challenges to the magnetar model. One example is the existence of the low magnetic field 
magnetar (Rea et al. 2010; Tong & Xu 2012). Alternatively, it is possible that magnetars are wind braking 
(Tong et al. 2012 and references therein). Wind braking would help us to explain the controversial timing 
results of the magnetar Swift J1822.3 — 1606. 

Swift J 1822. 3— 1606 is a magnetar candidate, discovered by Svw/r/BAT on 201 1 July 14 (Cummings 
et al. 2011). Up to now, different timing results are obtained for this source (Livingstone et al. 2011; 
Rea et al. 2012; Scholz et al. 2012). The reported period derivative differs by a factor about three. The 
corresponding characteristic magnetic field can be larger or smaller than the quantum critical field. This 
is directly related to whether this source is another low magnetic field magnetar or not. 

In Rea et al. (2012) and Scholz et al. (2012) observational papers, they mainly discussed the timing 
noise effect. In their opinion, it is the the large timing noise that results in different period derivative 
measurement of Swift J1822.3 — 1606. In this paper, we explore another effect. The period derivative of 
Swift J1822.3— 1606 may be decreasing with time. Therefore, it is natural that different period deriva- 
tives are obtained using different data sets. The physical reason may be that magnetars are wind braking 
(Tong et al. 2012). A decaying particle wind after outburst will result in a decreasing period derivative. 

Model description and quantitative calculations are presented in Section 2. Discussions and conclu- 
sions are presented in Section 3. 
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2 MODELING THE SPIN DOWN RATE OF SWIFT J1822.3 1606 

2.1 Description of observations and theory 

In Rea et al. (2012), they reported two period derivatives of Swift J1822.3 — 1606 (Section 3.2 there). 
Using the first 90 days observations, a period derivative P = 1.6(4) x 10~ 13 is obtained (the last digit 
uncertainties are at 1 a confidence level). Considering the whole 275 days data, the corresponding period 
derivative is P = 0.83(2) x 10 -13 . These two values provide some hints that the period derivative is 
decreasing with time. The large uncertainty in short time data set may be caused by timing noise. Similar 
behavior can also be seen in Livingstone et al. (2011) and Scholz et al. (2012). In Livingstone et al. 
(2011), using 84 days observations, a period derivative P = 2.55(22) x 10~ 13 is reported. Using 402 
days observations, Scholz et al. (2012) reported three solutions of period derivatives: P = 0.683(21) x 
10~ 13 (fitting with period and period derivative), P = 1.71(7) x 10~ 13 (fitting with period and two 
period derivatives), P = 3.06(21) x 10~ 13 (fitting with period and three period derivatives). 

Similar behaviors are also seen in other magnetars. Since the early stage of magnetar timing studies, 
it is found that magnetars have a higher level of timing noise than normal pulsars (Gavriil & Kaspi 
2002; Woods et al. 2002). Large period derivative variations are seen in AXP IE 2259+586 (Kaspi et al. 
2003), AXP IE 1048.1-5937 (Gavriil & Kaspi 2004), SGR 1806-20 (Woods et al. 2007), and AXP IE 
1547.0-5408 (Camilo et al. 2008). Two neat examples are AXP XTE J1810 197 (Camilo et al. 2007) 
and the radio loud magnetar PSR J1622— 4950 (Levin et al. 2012). In these two sources, a decreasing 
period derivative is observed while the star's X-ray luminosity is decreasing after outburst. Therfore, 
from previous observations, there may also be large timing noise in Swift J1822.3 — 1606. At same time, 
its period derivative may also decrease with time (a decreasing X-ray luminosity is also observed). This 
may explain why a lower period derivative is obtained when using longer time span of observations. 

The physics for a varying period derivative may be that magnetars are wind braking (Tong et al. 
2012). The decay of strong magnetic field will power the star's X-ray luminosity. At the same time, 
a (magnetism-powered) particle wind will also be generated. The rotational energy of magnetars will 
mainly be carried away by this particle wind. A varying particle wind will naturally result in a varying 
period derivative. The fluctuations of this particle wind may account for the large timing noise in mag- 
netars. Since both the X-ray luminosity and the particle wind luminosity are from magnetic field decay, 
a model independent estimate of the particle wind luminosity is L p ~ L x , where L p and L x are the 
particle wind luminosity and the X-ray luminosity, respectively. The origin for this particle wind may be 
either internal (e.g., low amplitude seismic activities, Thompson & Duancan 1996), or magnetospheric 
(e.g., coronal particles, Beloborodov & Thompson 2007). For details of wind braking of magnetars and 
discussion of other models, see Tong et al. (2012) and references therein. 

2.2 Calculations for Swift J1822.3 1606 

X-ray observations of Swift J1822.3 — 1606 have given its flux evolution with time. Using the flux evo- 
lution function and its extrapolations, we can calculate the theoretical period derivative as a function 
of time. The longest time span of X-ray observations of Swift J1822.3 — 1606 is done by Scholz et al. 
(2012, 400 days observations). According to Scholz et al. (2012), a double exponential flux decay model 
is prefered. 

F(t) = Fx exp[-t/n] + F 2 cxp[-t/r 2 ] + F q , (1) 

where F(t) is the 1 — 10 kcV source flux as a function of time, t is in units of days since BAT trigger 
time (MJD 55756.5), t\ = 15.5 days and t 2 = 177 days are the two decay time scales, F\ — 20.9 x 
10~ n ergcm~ 2 s _1 and F 2 = 1.74 x 10~ n ergcm~ 2 s _1 are the two flux normalizations, F q = 

3 x 10~ 3 x 10 -11 ergcm~ 2 s _1 is the fixed quiescent flux (constrained by ROSAT). See Scholz et al. 
(2012, Section 3.3 there) for details. 

1/2 

The rotational energy loss rate due to an isotropic particle wind is proportional to L p (Section 3 
in Tong et al. 2012). Therefore, the period derivative will evolve with time as (short term evolution, e.g. 
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several years) P(t) oc L^ 2 oc lU" 2 oc Fit) 1 / 2 . Including a constant factor 

P(t)=N F(t) 1 ' 2 , (2) 

where No is the normalization constant. The observational reported period derivative is the average 
value over a certain time span. Expanding the period at epoch t\, 

P(t)=P(t 1 ) + P(t 1 )(t-t 1 ), (3) 

where P(t) and P(ti) are the rotation period at time t and t\, respectively, P(ti) is the period derivative 
at t 1 . Therefore, the observational reported period derivative for time span t — 1\ is (t is the end time, t\ 
is the starting time) 

Pobs(t - it) = P(ti) = -^—(P(t) - P(h))- (4) 
t — ti 

Rewriting the above equation, 



Pob s (t-h) = — — / P(t')dt' (5) 

I - El Jti 

= N *T^T f F(t') 1/2 dt' (6) 
t-tx J tl 

= N g(t,h), (7) 



where g(t,ti) = J i 7, (i') 1 / 2 di'/(t — ti), F(t) is obtained by fitting the observational flux decay 
(equation (HJ). 

The timing of Livingstone et al. (2011) is done for time span 85.5 — 1.5 days since BAT trig- 
ger time. While the timing of Scholz et al. (2012) is for time span 404.5 — 2.5 days since BAT trig- 
ger time. According to equation (0, the ratio of period derivative between Scholz et al. (2012) and 
Livingstone et al. (2011) should be g(404.5, 2.5)/g(85.5, 1.5) = 0.48. The observational values are 
0.683(21)/2. 55(22), 1.71(7)/2.55(22), 3.06(21)/2.55(22), for solution 1, solution 2, and solution 3, 
respectively. The observation and theory are consistent within uncertainties. The same can also be done 
for the timing of Rea et al. (2012). Using the observational flux decay there, the theoretical period 
derivative ratio between 275 days and 90 days observation is 0.60. While the observational value is 
0.83(2)/1.6(4). The observation and theory are consistent with each other. 

We can also plot the theoretical period derivative as a function of time span. Employing the period 
derivative of solution 1 in Scholz et al. (2012) as the normalization, the predicted period derivative as a 
function of time is 

P ohs (t- ^ = ^(404.5-2.5)^1^. (8) 

The timing solutions in Livingstone et al. (2011), Rea et al. (2012), Scholz et al. (2012) are taken at 
different epoch (i.e., different ti). However, the differences are only one or two days. Therefore, this 
difference is negligible, ti = 2.5 is assumed in the following calculations (the ti value in Scholz et 
al. (2012)). Figure Q] shows the theoretical period derivative and the current observational data. The 
theoretical curve (using solution 1 in Scholz et al. (2012) as normalization) is consistent with the timing 
of Rea et al. (2012). The large uncertainties in the timing of Livingstone et al. (201 1) and 90 days timing 
result of Rea et al. (2012) may be due to timing noise. 

In the future, when longer time span observations are available, a smaller period derivative is ex- 
pected. For example, 800 days of timing observations will result in a period derivative P = 0.44 x 10" 13 . 
This is the theoretical period derivative averaged over 800 days. If separate timing can be done for the 
early 400 days and the late 400 days, a smaller period derivative is expected. Current 400 days tim- 
ing gives a period derivative P = 0.683 x 10~ 13 . A period derivative P = 0.19 x 10~ 13 is expected 
for the late 400 days timing only. It is about three times smaller. Future timing observations of Swift 
J1822.3 — 1606 will help us make clear whether its period derivative is decreasing with time or not. 
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Fig. 1 Theoretical period derivative as a function of observational time span. The continuous 
line is the theoretical period derivative. Red circles are timing data from Rea et al. (2012). 
Blue squares are timing data from Livingstone et al. (2011) and Scholz et al. (2012). The 
filled blue square is taken as normalization of the theoretical curve. The error bars are 3 a. 



3 DISCUSSIONS AND CONCLUSIONS 

The above calculations are mainly base on equation ([3j. In equation (0), only the first period derivative 
is included in the expansion. The observed P is the average value of period derivative over the ob- 
servational time span. During timing studies higher order period derivatives may also be included (e.g., 
solution 2 and solution 3 in Scholz et al. 2012). When higher order period derivatives are considered, the 
corresponding P will approach its instantaneous value at the expansion epoch. Therefore, the reported 
P represents earlier value when higher order period derivatives are included. If the physical spin down 
rate is decreasing with time, we should see a larger P when higher order period derivatives are included. 
This is just the three timing solutions in Scholz et al. (2012). Therefore, the three timing solutions of 
Scholz et al. (2012) provides us another evidence that the period derivative of Swift J1822.3 — 1606 is 
decreasing with time. 

When calculating the theoretical spin down rate, the particle wind luminosity is assumed to equal to 
the soft X-ray luminosity. The actual wind luminosity may have a more or less different value. After the 
outburst, the star's X-ray luminosity decreases with time. Since the particle wind is also from magnetic 
field decay, then it is natural that the wind luminosity also decreases with time. Therefore, a decreasing 
period derivative is always expected irrespective of the details of particle wind luminosity. In the long 
term run, the X-ray luminosity will return to its quiescent value. The particle wind will also relax to 
its quiescent state. The long term predicted period derivative is very sensitive to the condition of the 

1 12 

quiescent state. When assuming L p = L K , the period derivative at late time will be P oc F q , where 
F q is the quiescent flux. A quiescent flux ten times higher, the late time period derivative will be be three 
times larger. 

The surface dipole field obtained by assuming magnetic dipole braking is only the effective field 
strength. In the presence of strong particle wind, the rotational energy loss rate is amplified. For a given 
period derivative, the resulting dipole field will be much lower (Tong et al. 2012). In the actual case, the 
geometry (e.g., the magnetic inclination angle) will also affect the spin down history of the neutron star 
(Tong & Xu 2012). In the case of normal pulsars, the magnetic dipole braking assumption is a not too 
bad lowest order approximation (Xu & Qiao 2001). However, in the case of magnetars, the magnetic 
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dipole braking assumption will be too simple even to the lowest order approximation. An alternative 
is that magnetars are wind braking (Tong et al. 2012; Tong & Xu 2012). A decaying particle wind can 
result in a decreasing period derivative of Swift J1822.3 — 1606. 

Another explanation for the decreasing period derivative is the twisted magnetosphere model 
(Thompson et al. 2002; Beloborodov 2009). After the outburst, the magnetar magnetosphere gradu- 
ally untwists. Therefore, the effect dipole magnetic field will decrease. This will cause a decreasing 
period derivative. However, the twisted magnetosphere model may have difficulties in explaining the 
short time scale period derivative variations (Camilo et al. 2007; Levin et a. 2012). In the above wind 
braking of magnetars, the wind luminosity can vary dramatically on short time scales. Such difficulties 
no longer exists in the wind braking model. 

In conclusion, the different timing results of Swift J1822.3 — 1606 are caused not only by its timing 
noise, but also by its decreasing period derivative. The decreasing period derivative and large time noise 
may be both originated from wind braking. Future timing observations of Swift J 1822. 3 — 1606 will help 
to make clear whether its period derivative is decreasing with time or not. This would also help us to 
answer whether wind braking is important in this source or not. 
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